Spherical to deformed shape transition in the nucleon-pair shell model 
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A study of the shape transition from spherical to axially deformed nuclei in the even Ce isotopes 
using the nucleon-pair approximation of the shell model is reported. As long as the structure of 
the dominant collective pairs is determined using a microscopic framework appropriate to deformed 
nuclei, the model is able to produce a shape transition. However, the resulting transition is too 
rapid, with nuclei that should be transitional being fairly well deformed, perhaps reflecting the need 
to maintain several pairs with each angular momentum. 
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I. INTRODUCTION 

One of the most fascinating aspects of nuclear struc- 
ture is the existence of competing modes of collective 
motion, with phase transitions linking them. This has 
been discussed extensively within the context of the in- 
teracting boson model (IBM) |l| , with the vertices of the 
Casten triangle representing extreme modes of collective 
motion in the ground state 0. 

There has been much effort to develop theories that are 
able to describe the competing modes of nuclear collec- 
tivity and the associated transitions from one to another. 
This is a major theme of the ongoing effort to develop 
a unified density functional theory of atomic nuclei. It 
has also been a major theme of the phenomenological 
IBM, which has recently seen a significant microscopic 
advance through the development of a mapping proce- 
dure that translates the physics of density functionals to 
the IBM parameters in the various collective regions 

A model that is closely linked in spirit to the IBM is the 
nucleon-pair approximation (NPA) of the shell model (of- 
ten called the nucleon-pair shell model in the literature) 
0,0]. Like the IBM, the NPA is based on a truncation in 
terms of pair degrees of freedom. Unlike the IBM, how- 
ever, the microscopic structure of the pairs is preserved 
throughout. 

Recently, it has been shown that the NPA truncated 
to the dominant S and D pair degrees of freedom is able 
to qualitatively reproduce all of the phase transitions of 
the IBM, including the transition from vibrational U(5) 
behavior to axially symmetric SU(3) rotational behavior 
These studies were of a model nature, however, not 
making contact with real nuclear systems. 

In this work, we return to the issue of whether the NPA 



is able to describe the shape transition from vibrational 
to rotational behavior in real nuclei. We carry out the 
analysis for the even-^4 Ce isotopes from 142 Ce through 
148 Ce, where such a shape transition has been observed 
experimentally. We find that when we use the traditional 
approach for generating the dominant collective pairs of 
the model, we cannot reproduce the shape transition. We 
then suggest the use of an alternative approach to the 
structure of these pairs for rotational nuclei that incor- 
porates more of the physics that builds their correlations. 
With this new approach to the structure of the collective 
pairs in the rotational region, we are able to achieve the 
observed shape transition. 

The structure of this article is as follows. In Sec. UTJ we 
briefly review the basic ingredients of the NPA, including 
a discussion of the traditional approach for defining the 
structure of the collective pairs to be used when applying 
the method. In Sec. IIII| we describe the alternative 
formalism for building the collective pairs appropriate to 
deformed nuclei. In Sec. IIV1 we discuss other input to 
our calculations for the even-mass Ce isotopes, for which 
the results are presented in Sec. [V] Finally, in Sec. I VII we 
summarize the key conclusions of the work and outline 
some issues for future consideration. 



II. THE NUCLEON-PAIR APPROXIMATION 
OF THE SHELL MODEL 

The nucleon-pair approximation of the shell model 
(NPA) is a truncation strategy for the nuclear shell model 
based on the dominance of a few selected coherent pair 
configurations. When viewed as a microscopic realization 
of the IBM, the selected nucleon pairs are the lowest spin- 
zero and spin- two pairs for neutrons and protons, called S 
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and D pairs, respectively. The NPA does not map these 
collective pairs onto bosons, however, working in terms of 
the collective pair configurations directly. Furthermore, 
the NPA does not restrict the collective pairs that are 
included to the lowest spin-zero and spin-two pair con- 
figurations for neutrons and protons only, permitting in 
principle other pair configurations to be included. 

The NPA traces back to the generalized seniority 
scheme [3] and the broken pair approximation (BPA) Q. 
These methods could not be used as a practical and gen- 
eral truncation strategy for the shell model, however, un- 
til the development by Chen 0] in 1997 of a set of recur- 
sive formulae that enabled the calculation of the matrix 
elements of realistic nuclear hamiltonian in a basis built 
up in terms of these collective pairs. We italicize the term 
realistic because of the computational limitation to date 
of the methodology to hamiltonians involving multipolc 
pairing interactions in the like-nucleon channels and sep- 
arable multipole-multipole interactions in the neutron- 
proton channel. With such assumed semi-realistic hamil- 
tonians and codes based on the Chen's recursion formu- 
lae, it is now possible to treat systems involving several 
pairs of each type of particle and several collective pair 
degrees of freedom. 

There are several key steps in the NPA strategy. One 
involves a choice of the effective shell-model hamiltonian, 
along the lines just discussed. The second involves an as- 
sumption of the pair degrees of freedom to be included, as 
dictated by the physics of the problem under discussion 
and computational limitations. Another key ingredient 
is the assumed structure of the collective pairs that are 
retained. The standard approach has been to obtain the 
lowest S pair from a number-conserving BCS treatment, 
or equivalently a variational description in terms of states 
with generalized seniority zero, for neutrons and protons, 
respectively. The prescription typically used for the low- 
est excited pairs is a variational treatment of these pairs 
in terms of generalized seniority- two states, using the S 
pairs already determined. An alternative approach some- 
times used for the structure of the D pair is to generate 
it from the commutator of the quadrupole operator with 
the S pair. The traditional prescription is reasonable for 
systems that are meaningfully described by the general- 
ized seniority or broken pair scheme, namely for cither 
semi-magic or near semi-magic nuclei. Once there are a 
large enough number of both types of particles, however, 
the expectation is that generalized seniority ceases to pro- 
vide a useful description of the states of the system. In 
such cases, an alternative to the traditional approach for 
choosing the dominant collective pairs must be sought. 
Indeed, such considerations have led to the difficulties 
that have persisted for many years in microscopically de- 
riving the IBM for all but spherical C/(5)-like nuclei, at 



least until the recent work in which the collective surface 
emerging from density functional treatments is mapped 
directly into a boson space, with no direct consideration 
of the internal structure of the collective pairs 0] . 

In the next section we discuss how to build more ap- 
propriately the dominant collective pairs for our study of 
the shape transition from vibrational to deformed nuclei 
using the NPA, and then in the results that follow in Sec- 
tion IV we show that the choice of this new strategy for 
building the collective pairs, a strategy borrowed from an 
earlier effort to build the relevant collective pairs in de- 
formed nuclei, is critical to obtain a meaningful descrip- 
tion of the shape transition from vibrational to deformed 
nuclei in the region we study, namely the even-mass Ce 
isotopes. 

III. MICROSCOPIC STRUCTURE OF 
COLLECTIVE NUCLEON PAIRS 

As we will see in next section, it is critical to use collec- 
tive pairs that are tuned to the dynamics of the collective 
region under investigation. In the vibrational region, the 
traditional pairs that have been used in the NPA are most 
likely appropriate as they follow from the dynamics of 
generalized seniority. In deformed nuclei, however, such 
pairs do not incorporate any effects of the proton-neutron 
interaction. Thus, we propose to use a prescription first 
discussed by Pittel and Dukelsky [9] in the context of the 
microscopic IBM for deformed nuclei. Their prescription 
derives from the use of Hartree-Fock-Bogolyubov (HFB) 
approximation to variationally generate an optimum in- 
trinsic pair, from which pairs of definite angular momen- 
tum can then be projected. Similar considerations to ob- 
tain coherent SD nucleon pairs appropriate to the NPA 
calculations based on the field approximation were stud- 
ied in Ref. [13 ■ 

We begin the procedure by carrying out an axially- 
symmetric HFB calculation in the same space and with 
the same hamiltonian that is used in the NPA study. 
Once the HFB calculations have been carried out, we 
then transform to the canonical basis, namely the basis 
in which the density matrix p*-™ 1 ' (r) for particles of type 
t [y or 7r) and projection m, is diagonal. Letting Tj m 
denote the operator that creates a nucleon of type r in 
orbit (jm), we can express the single-particle creation 
operators in the canonical basis by 

T am = (t) Tj m , (1) 

3 

where (r) is the matrix that transforms the density 
matrix to diagonal form for projection m and particle 
type r. The intrinsic HFB wave function in the canonical 
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basis can be expressed as 
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is a coherent pair creation operator and |0) denotes the 
doubly-magic core. 

The quantities entering ([3]) can be obtained directly 
from the diagonal elements of the density matrix, viz: 
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The correlated pair creation operators T\ contain com- 
ponents with all allowed even angular momenta, 



n = E4 r) 4 L)t 



(4) 



The square of the amplitude a L determines the impor- 
tance of the correlated pair of type r with angular mo- 
mentum L in the intrinsic wave function and is given by 
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The structure of the correlated pair If with angular 
momentum L is given by 
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and is the information of interest for use in the NPA when 
applied to well-deformed nuclei. 



TABLE I: Single-particle energies (in MeV) used in the cal- 
culations for the Ce isotopes described in the text. They are 
based on low- lying energy levels of 133 Sb and 133 Sn. 





Sl/2 


^3/2 


^5/2 


97/2 


hn/2 






2.990 


2.708 


0.962 


0.000 


2.793 






Pl/2 


P3/2 


75/2 


fr/2 


hg/2 


*13/2 




1.656 


0.854 


2.005 


0.000 


1.561 


1.800 



remaining protons over the orbits of the 50 — 82 major 
shell and the remaining neutrons over the orbits of the 
82 — 126 major shell. 

The hamiltonian, as dictated by the computational 
structure of the NPA code, takes the form 



h= E E'-:- E ^4 (L) -4 J 

T = f,TT \ j L = 0,2 



(5) 



It includes a single-particle energy term and a two-body 
interaction consisting of multipole pairing interactions 
in the nn and pp channels and a separable quadrupole- 
quadrupole interaction in the pn channel. 

The single-particle energies are extracted from the 
spectra of 133 Sb and 133 Sn, and are listed in Table HI The 
two-proton interaction parameters are obtained from the 
binding energy and 2f energy of 134 Te. The values that 
emerge are 



-0.18 MeV, G n2 = . 



Note that a pure monopole pairing interaction reproduces 
both the binding energy of 134 Te and the excitation en- 
ergy of the 2^ state. The two-neutron parameters are ob- 
tained from the binding energy and 2^~ energy of 134 Sn. 
The resulting values are 

G„ = -0.13 MeV, G v2 = -0.012 MeV . 

The strength of the proton-neutron quadrupole interac- 
tion is obtained by requiring an optimum description of 
both the binding energy and the lowest 2 + excitation en- 
ergy of 136 Te, which has two valence protons and two 
valence neutrons, and is given by 

K = -0.20 MeV . 



IV. PARAMETERS IN OUR NPA 
CALCULATIONS 

We treat the even-mass Ce isotopes by assuming a 
Z = 50, N = 82 doubly-magic core and distributing the 



V. CALCULATED RESULTS 

In this section, we present our calculated results ob- 
tained using the NPA for the even-mass Ce isotopes from 
142 Ce through 148 Ce. That this series of isotopes exhibits 
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TABLE II: Ratios of experimental excitation energies of the 
Ce isotopes from Ref. [111 ]. R4/2 gives the ratio of the 4^ to 
the 2+ excitation energies, and Rq/2 the ratio of the 6^ to 2f 
energies. The last column, denoted SU(3), gives the results 
in the rotational limit. 
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FIG. 1: Deformation energies (iJdcformcd — -Esphcricai) from 
axially-symmetric HFB calculations for the even-mass Ce iso- 
topes as a function of the number of valence neutrons N v 



a shape transition from vibrational to rotational nuclei 
can be seen from Table [TT1 where we present the ratios 
of excitation energies found experimentally for the 4^~ 
and 6f levels relative to the 2f levels through this chain. 
The nucleus 142 Ce, with N v = 2 valence neutrons, has 
a vibrational character, with a R4/2 ratio of 1.90, while 
148 Ce with N v — 8 valence neutrons has a fairly rota- 
tional character with R4/2 = 2.87 and Rq/2 = 5.32. The 
transition from vibrational to rotational is fairly smooth. 

In Fig. [TJ we show the deformation energies 
(^deformed - Spherical) that derive from our HFB cal- 
culations of these nuclei. As expected, the deformation 
energy grows gradually with increasing neutron number, 
with the gain in energy reaching a value of 1.27 MeV in 
148 Ce. This is not as large a deformation energy as we 
would have liked to have seen for a well-deformed nucleus, 
perhaps reflecting the fact that our model space does not 
include all the orbits that could contribute meaningfully. 
Nevertheless, the results are suggestive that the HFB cal- 
culations are producing the observed shape transition. 

In Fig. [2j we present the calculated HFB results for 
|cjl| 2 for these isotopes, illustrating the gradually increas- 
ing importance of L 7^ pairs in the ground band as the 
deformation grows. By the time we reach 148 Ce, D pairs 
contribute to the intrinsic state of the ground band as 



FIG. 2: Contributions of pairs with definite angular momen- 
tum to the HFB intrinsic states calculated for the Ce isotopes 
with different neutron pair number N v . 



much as S pairs. 

Next we present the results of our calculations using 
the NPA. In all calculations our model space includes the 
lowest S pair, the lowest D pair and the lowest G pair 
(with spin four), both for neutrons and protons. How- 
ever, because of computational considerations we only 
permit a single G pair for each type of particle. Fig. 
[3] shows several sets of calculated results in comparison 
with the experimental data. The spectra denoted "BPA" 
are based on the use of the traditional broken pair ap- 
proximation (BPA) prescription for the S, D and G pairs, 
as described in Section II. The spectra denoted "HFB" 
are based on the use of pairs obtained from axially- 
symmetric HFB calculations, as described in Sec. IIIII 
The spectra denoted "COM" include the S and D pairs 
from the two prescriptions in the same calculation, with 
their nonorthogonality taken into account. We return 
later to our reason for including the latter results. Re- 
sults are shown up to 5.5 MeV in excitation energy, both 
for the experimental and calculated spectra. 

There are several clear messages that emerge from 
these results. First, the traditional broken pair approx- 
imation (BPA) results are unable to describe the shape 
transition from vibrational to rotational as we proceed 
through the Ce chain. This is not surprising, consider- 
ing our earlier motivating remarks. The HFB results, in 
contrast, do give rise to a clear shape transition. Though 
142 Ce does not exhibit a rotational pattern, the addition 
of just two neutrons leads to a fairly robust rotational 
spectrum, which becomes even more pronounced with 
the addition of further neutrons. 

These conclusions are made even more striking when 
we show the energies of the associated gamma rays 
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FIG. 3: Comparison of the spectra calculated using the NPA for the even Ce isotopes (panels b-d) with the corresponding 
experimental spectra (panel a, denoted by "EXP" ) . The model spaces are constructed using generalized seniority-two states for 
both protons and neutrons (denoted by "BPA"), and SDG- nucleon pairs obtained from HFB calculations (denoted by "HFB"), 
in panels (b) and (c), respectively. The model space of panel (d) is a direct sum of the two bases in panels (b) and (c), and is 
denoted "COM". See the text for details. 



emerging from the two sets of calculations, which are 
given in Fig. @] The HFB results show a rapid emer- 
gence of a deformed pattern (roughly parallel lines) with 
increasing neutron number, whereas none is seen in the 
BPA spectra. 

Interestingly, the BPA results do not even give a good 
description of 144 Ce, with only four valence neutrons. 
The excitation energy it predicts for the lowest 2 + state 
is significantly larger than that seen experimentally. This 



is to be contrasted with the HFB results, which already 
gets a fairly good description of the lowest 2 + excitation 
energy in 144 Ce. 

While we have focussed on energy spectra in Figs. [3] 
and HJ and in the discussion that followed, we can get a 
clearer picture of why it is important to use dynamically- 
determined collective pairs by looking at the absolute 
binding energies of the ground states that emerge in these 
calculations with respect to the 132 Sn core. These results 
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FIG. 4: Comparison of associated gamma rays [E(I) — E(I — 
2)] emerging from the three sets of calculations with the cor- 
responding experimental data. Panels (a-d) have the same 
meaning as in Fig. 3. 
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FIG. 5: Binding energies arising from the NPA calculations 
for the even Ce isotopes as described in the text. The curve 
denoted BPA refers to the use of the broken pair approxi- 
mation to define the collective pairs; the curve denoted HFB 
refers to the use of the Hartree-Fock-Bogolyubov approxima- 
tion to generate them. 



are shown in Fig. [5j When we use the broken pair ap- 
proximation prescription for the collective pairs of the 
model we obtain much less binding than with the HFB 
pairs, which include effects of the proton-neutron inter- 
action in their structure. While this is especially true 
for the most well-deformed nuclei, it is also true, albeit 
to a much less extent, in 142 Ce, where the use of HFB 
collective pairs produces roughly 0.1 MeV in additional 
binding. 

While the method is able to produce a shape transition 
from vibrational to rotational nuclei, it produces it much 



too rapidly. As is evident both from Fig. [3]and even more 
so from Fig. 21 the model already produces a fairly well- 
deformed spectrum for 144 Ce, which according to Table 
HTl is clearly transitional. 

It is not too surprising that we are unable to ade- 
quately describe transitional nuclei, since our axially- 
symmetric HFB solution is not very meaningful for those 
nuclei. The fairly small HFB deformation energy for 
those nuclei, as seen in the Fig. [TJ suggests that the HFB 
solution is not stable against zero-point fluctuations. One 
possible thought is that we might include more than one 
pair for each angular momentum in this region. As a 
first step in this direction, we have carried out calcula- 
tions in which we include both the HFB pairs and the 
BPA pairs, taking into account their nonorthogonality, 
to see whether this improves our description of the tran- 
sitional region. Those are the results denoted "COM" 
in both Figs. [3] and SJ The inclusion of a second S and 
D pair for each type of nucleon improves the description 
substantially, giving a better reproduction of the proper- 
ties of all nuclei considered. It markedly improves upon 
the HFB results for 142 Ce, placing the 2^ state much 
closer to its experimental location. However, as is evi- 
dent from Fig. |H the shape transition is still too rapid. 
Nevertheless, we believe that the idea of including more 
than one pair with each angular momentum is a possible 
approach for future consideration, as it can be readily 
included in the analysis, but further thought as to which 
pairs to include, especially in the transitional region, is 
needed. 

It is also interesting to note that even in the most de- 
formed nucleus considered, 148 Ce, the moment of iner- 
tia that derives from the best NPA calculations, those 
denoted "COM", is still too small and the spectrum is 
accordingly too spread-out. Interestingly the same result 
emerged in the first efforts to microscopically derive the 
IBM in deformed nuclei by mapping the fermionic density 
functional onto a corresponding one in the boson space. 
Recently it was shown that it is necessary to add a fur- 
ther term to the IBM hamiltonian of the form aJ(J + 1) 
to better reproduce the moment of inertia of the rota- 
tional systems that were studied in Ref. [lij . The NPA 
seems to also have need for an additional contribution 
to the rotational moment of inertia, which perhaps has a 
similar origin. It is not clear at this time, however, how 
to incorporate it in the NPA formalism. 



VI. SUMMARY AND CLOSING REMARKS 

In this paper, we have addressed the issue of whether 
or not the NPA is able to describe the transition from vi- 
brational to rotational nuclei in real nuclear systems. We 
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focus on the even-mass Ce isotopes, in which such a shape 
transition is known to occur experimentally. We see that 
with the use of the traditional collective pairs that have 
been used in previous applications of the model, namely 
pairs derived from the broken pair or generalized senior- 
ity approaches, the model is not able to describe such a 
shape transition. Such pairs leave out important effects 
of the proton-neutron interaction and thus are not suit- 
able in a truncated description of collective nuclei with 
active neutrons and protons. In contrast, a description 
that uses collective pairs from the HFB approximation 
does contain the requisite correlation effects and is able 
to produce the observed shape transition. On the other 
hand, the shape transition produced is too sharp, reflect- 
ing the fact that the correlated pairs produced in the 
HFB calculations are not adequate for the transitional 
region. Furthermore, the moments of inertia are not very 
well reproduced. 

The difficulty in obtaining appropriate collective pairs 
to use in the NPA in the various nuclear regions is rem- 
iniscent of the same difficulty that arose in early efforts 
to microscopically derive the IBM for these different re- 
gions. It is only with the recent development of methods 
that map the collective surface directly onto an associ- 
ated IBM collective surface that it has been possible to 
consistently obtain a microscopic derivation of the IBM. 



In the absence of an analogous procedure for the NPA, 
we have suggested the possibility of using several collec- 
tive pairs in those regions in which a single HFB pair for 
each angular momentum is not sufficient. Further work 
to identify how to optimally obtain those pairs in transi- 
tional regions is still needed. Likewise, we have seen that 
the moment of inertia that results from our best calcu- 
lation for well deformed nuclei is too small and further 
work is needed to find a way to increase the rotational 
moment of inertia in the NPA framework. 
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